Quasi-two-dimensional liquid ^He 
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Quantum Monte Carlo simulations at zero temperature of an ensemble of He atoms in the 
presence of smooth Mg, Li, Na and K substrates yield strong evidence of a thermodynamically stable 
liquid ^He monolayer for Na and K. The effective 2D density is 6' ~ 0.02 A~^ on Na, making it the 
lowest density liquid in Nature. Its existence is underlain by zero-point atomic motion perpendicular 
to the substrate, whose effect is that of softening the short-range repulsion of the helium interatomic 
potential. The monolayer film is expected to turn superfluid at a temperature Tc ~ 1 mK. No liquid 
film is predicted to form on Mg, or on stronger substrates such as graphite. 



PACS numbers: 67.30.ej,67.30.hr 

Decades of experimental and theoretical investigation 
have yielded evidence that the behaviour of sufficiently 
thin films of "^He closely mimics that theoretically pre- 
dicted for a system strictly confined to two spatial di- 
mensions (2D), atomic zero-point motion in the direction 
perpendicular to the subtstrate notwithstanding. For ex- 
ample, the superfluid transition in thin liquid ^He films 
has been observed [IHl] to conform to the 2D Kosterlitz- 
Thouless paradigm [5]. The effect of substrate corruga- 
tion and/or roughness is generally limited to the appear- 
ance of insulating phases, either disordered or crystalline, 
registered with the underlying substrate, at commensu- 
rate coverages. Thus, thin films of *He constitute an 
ideal quasi-2D system, on which fundamental properties 
of strongly interacting Bose fluids in two dimensions can 
be investigated experimentally. 

The lighter isotope of helium (^He) could in princi- 
ple play the same role for Fermi systems, as ^He atoms 
are composite spin-1/2 particles. However, while '^He 
forms registered solid monolayers on corrugated sub- 
strates, with atoms sitting at preferential adsorption 
sites, the conventional wisdom is that on any substrate 
strong enough to be wetted by it, '^He will form only 
a "thick" film, whose physical properties are essentially 
those of the bulk liquid phase. The theoretical argu- 
ment is that, unlike "^He, whose 2D equilibrium phase 
at temperature T—0 is a liquid with a binding energy of 
approximately 1 K per atom [3], no self-bound liquid is 
deemed to exist in 2D for ^He [7HTT]. 

A significant amount of experimental work has fo- 
cused on a possible quasi-2D liquid phase of highly dilute 
■^He, floating atop a thin "'He film adsorbed on differ- 
ent substrates P^^HTl] . While some of these studies have 
yielded some evidence of it, consensus is still lacking as to 
whether the results unambiguously point to the existence 
of a self-bound quasi-2D ■^He liquid [15]. Furthermore, 
interesting as the physics of a film of a mixture of the 
two isotopes of helium undoubtedly is [TB] , the physical 
environment experienced by a ^He atom moving on top 



of (or through) a superfluid ^He film cannot be regarded 
as fully equivalent to that of a planar surface, as the 
•^He atom can scatter off surface excitations (ripplons) of 
the underlying "^He film. Thus, the question remains of 
whether a substrate exists, on which ^He by itself could 
condense into a quasi-2D liquid. 

Recent experiments [T7] suggest that this may be the 
case for ^He adsorbed directly on graphite, specifically 
the apparent formation of liquid puddles in an adsorbed 
monolayer, at a coverage ^ 0.01 A~^, i.e., far less than 
that corresponding to a registered solid monolayer. In 
this case too, however, an unambiguous interpretation of 
the experimental evidence is complicated by the expected 
significant role played by surface imperfection and/or in- 
homogeneities [15] . 

In this Letter, we provide theoretical evidence that '^He 
will form quasi-2D liquid monolayers at T=0 on Na and 
K substrates, with a coverage of approximately 0.03 A~^ 
on K, and 0.02 A~^ on Na. We base this statement on 
atomic energetics computed by first principles Quantum 
Monte Carlo (QMC) simulations of a realistic model of a 
•^He adsorbate. We have considered four weakly attractive 
substrates, namely Mg, Li, Na and K (listed in order of 
decreasing attractiveness) . Our findings are that a liquid 
monolayer will only form on substrates that are suffi- 
ciently weak; our calculations cannot exclude the possi- 
bility of a very dilute liquid on Li, which appears to be the 
borderline case, whereas on the stronger Mg substrate, no 
thin '^He ffim will form. In all predicted self-bound liquid 
phases, zero-point motion in the direction perpendicular 
to the substrate is significant, and has the effect of soft- 
ening the hard-core repulsion of the interatomic '^He pair 
potential at short distances, in turn allowing for a self- 
bound liquid to occur. 

Our system of interest is modeled as an ensemble of N 
•^He atoms, {N/2 of either spin component, i.e., the sys- 
tem is unpolarized) regarded as point particles, moving 
in the presence of an infinite, smooth planar substrate 
(positioned at z = 0). The system is enclosed in a ves- 



sel shaped as a parallelepiped, with periodic boundary 
conditions in all directions. The nominal '^He coverage 
is 6 ^ N/A, A being the area of the substrate. The 
quantum-mechanical many-body Hamiltonian is the fol- 
lowing: 
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where m is the '^He atomic mass, V is the potential de- 
scribing the interaction between two helium atoms, only 
depending on their relative distance, and U is the poten- 
tial describing the interaction of a ^He atom with the sub- 
strate, also depending only on the distance of the atom 
from the substrate. We use the accepted Aziz poten- 
tial [IH] to describe the interaction of two ^He atoms. 
For each of the substrates studied here, the U term in 
Eq. (II]) is a potential proposed by Chizmeshya, Cole 
and Zaremba [19j. The assumption of a smooth, planar 
substrate is clearly an important one; its justification is 
provided by the relative weakness of the substrates con- 
sidered here. 

We investigate the ground state of (IT]) by means of 
a QMC technique, specifically Diffusion Monte Carlo 
(DMC), which projects the lowest-energy component out 
of an initial trial wave function ^y. This is a well- 
established methodology HHj: which has been utilized to 
study ground state properties of a wide variety of quan- 
tum many-body systems. The trial wave function that 
we utilized in this work has the form ^P — _F$, where, 
using standard notation (see, for instance, Refs. pTll^ l 



exp 



tN 



N 

E^ 



N 

E 

i<j<k=l 



Wi-jk) 



(2) 



and $ = {n:Li /(z.)} ^{t)Hi), 0(a) (" =1,1) given by 



Det, 



exp 



iki 



N 



kjj^kj 



(3) 



where s^ = (xi,j/i), the wave vectors k^ span the N/2 
allowed momenta of a 2D Fermi sea; in Eq. (JS]), i and j 
range from 1 to N/2 {N/2 + 1 to N) for a =t (i). The 
two-, three-body and backflow correlation functions u, w 
and ry, as well as the single-particle orbital f{z) are opti- 
mized based on the procedure described in Ref. [23] . 
We make use of the well-known "fixed-node" approxima- 
tion (FNA) in order to circumvent the inevitable sign 
problem that affects any fermion QMC method. We are 
mainly interested in the energy per '^He atom computed 
as a function of coverage, and the FNA has been shown 
to provide accurate, variational energy estimates for bulk 
3D ^He [23]. The simulations carried out in this work, 
for which results are shown, are for a system of A^=26 
•^He atoms; the DMC time step is 10~^ K~^, and the 
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FIG. 1. (Color online) Ground state energy per atom e(Q) (in 
K) computed by DMC for a ^He film adsorbed on different 
substrates, as a function of coverage Q (in A~^). The bind- 
ing energy of a single atom on the given substrate is taken as 
the reference (zero) value for e(ff). Open symbols represent a 
second set of data for the Li substrate, obtained by applying 
estimates for the finite-size, time step, population control and 
fixed-node approximation errors. Solid lines are spline inter- 
polations. The data labelled "2D" refer to strictly 2D '^He. 
Statistical errors are smaller than the size of the symbols. 



population comprises 2,000 random walkers. We discuss 
systematic errors affecting our calculation after illustrat- 
ing our results. 

Figure IT] shows the ground state energy per atom e{Q)^ 
computed by DMC as a function of coverage on the dif- 
ferent substrates considered here. In all cases, the refer- 
ence energy is taken to be the binding energy of a single 
•^He atom on the given substrate [52]. Also shown for 
comparison is the result for the strictly 2D system. The 
existence of a thermodynamically stable quasi-2D liquid 
phase is signalled by a local minimum of e(0). One can 
clearly see that, starting from the strictly 2D case, and 
as the substrate is rendered weaker, i.e., atomic motion 
in the direction perpendicular to the substrate becomes 
more significant, the curve e(0) bends downward. On Mg, 
namely the most attractive of the four substrates, eiff) re- 
mains monotonically increasing, i.e., no stable liquid film 
forms. The corresponding curve for any substrate more 
attractive than Mg (e.g., graphite), lies between that for 
Mg and the one labelled 2D in Figure [T] implying no 
liquid "^He film on graphite as well, in conflict with the 
claim of Ref. [IT], under the assumption of a smooth 
substrate. 

The downward bending of e(ff) becomes increasingly 
pronounced on Li, Na and K, and a clear minimum can 
be observed on the last two substrates. On K, it is lo- 
cated at 6'=:0.03 A~^, whereas on Na at a coverage 0=0.02 
A^^. On Li, which is the weakest substrate on which "^He 
will form a superfluid monolayer |26| , our numerical data 
show no evidence of a stable liquid '^He film, although. 



as discussed below, Li may be a borderline case, and we 
cannot exclude a weakly self-bound liquid at a coverage 
close to 0.015 A-2. 




FIG. 2. (Color online) ''He density profile n{z) (arbitrary 
units) in the direction perpendicular to the surface, on a Mg, 
Li, Na and K substrate. The height of the main peak is greater 
the stronger the substrate (see text). Results shown are for a 
coverage 0=0.03 A~^. 



In order to gain insight into the physical properties of 
the adsorbed liquid films, we compute different structural 
properties, using extrapolated DMC estimators, typically 
adopted in DMC to obtain ground state expectation val- 
ues of quantities other than the energy |57]. Although 
they are not numerically exact, as they retain some of 
the bias of the starting trial wave function, these estima- 
tors are usually fairly reliable, especially if, as observed 
in this work, the difference between the variational and 
"mixed" estimate is small. 

Fig. p] shows the integrated "^He density profile n{z) = 
J dxdy p{x,y,z), where p{x,y,z) is the 3D ^He density, 
for the four substrates considered here. The height of 
the main peak is greater the stronger the substrate. Re- 
sults shown are for a coverage 9 = 0.03 A~^, but we 
have found this quantity to be nearly independent of 9, 
i.e., the density profile is almost entirely determined by 
the interaction of the helium atoms with the substrate, 
not surprisingly given the relatively large average planar 
interatomic separation in the range of coverage consid- 
ered here. The difference in the spread of the '^He atomic 
wave function in the z direction between the Mg sub- 
strate and the other ones, as well as the greater distance 
from the substrate, is evident. In particular, for the two 
substrates for which a self-bound quasi-2D liquid exists, 
the rms average excursion of each atom from its most 
probable distance from the substrate is of the order of a 
third of the average planar interparticle distance. 

With such a large spread, one might wonder to what 
extent the 2D characterization might be appropriate for 
such a system. A useful quantity is the angularly av- 
eraged, "reduced" pair correlation function g{R), with 







FIG. 3. (Color online) Reduced pair correlation function g{R) 
(Eq. p])), for ^He in 2D (sohd line), on a Mg substrate and 
on a K substrate, for a coverage d = 0.03 A~'^. The curve for 
which g{0) is finite is for the K substrate. 



R — y'x^ + y^ and 

5(^' y) = -J(p / c?2;'d2/' n{x + x',y + y') n{x' , y') (4) 

with n(a;, y) — J dz p{x, y, z). The more two-dimensional 
an adsorbed film, the more closely g{R) mimics the pair 
correlation function of a strictly two-dimensional system 
of the same coverage. Fig. 3 shows results for g{R) com- 
puted in the strictly 2D case, on a Mg substrate and on 
a K substrate, at a coverage 9 = 0.03 A~^. The results 
for the 2D case and on Mg are very similar, the main 
peak being slightly higher in 2D, and ^He atoms being 
able to come to slightly closer planar distances on a Mg 
substrate. In other words the physics on a Mg substrate 
approaches fairly closely the 2D limit. On the other hand, 
the result on a K substrate is qualitatively different, the 
main feature being a finite value of g{R) in the i? — > 
limit. This is a direct consequence of the large zero-point 
motion in the direction perpendicular to the substrate, 
whereby ^He atoms can occasionally be "on top" on each 
other. For this reason, the qualified "quasi" seems par- 
ticularly appropriate when referring to this liquid phase 
as 2D. 

Thus, as also suggested in Ref. [55], the main physical 
result of the large zero-point excursion is that of effec- 
tively softening the short-distance repulsive core of the 
interaction between two ^He atoms. By rendering the 
pair potential less repulsive, the substrate allows for the 
formation of a self-bound liquid. Based on this argu- 
ment, we predict that self-bound liquid phases of '^He 
should exist at T = on Rb and Cs substrates as well, 
at higher coverages than on K (the binding energy of a 
single ^He on Cs and Rb substrates is numerically very 
close to that on K) . This confirms an early prediction by 
Carraro and Cole [33] that ^He can wet substrates that 
are not wetted by ''He, in this case Cs. Conversely, and 
in some respects counter-intuitively, on a stronger, more 



attractive substrate such as Mg, or graphite, the physics 
of the system closely approaches the 2D limit, in which 
no self-bound liquid is expected to exist. 

We now discuss the main sources of systematic error 
of our calculation, in order to assess their expected effect 
on the physical conclusions outlined above. These are, 
for the energy, the finite size of the simulated system, 
the DMC time step error, the bias due to a finite walker 
population and the FNA, based on the nodal structure 
given by Eq. M |30j . An important general remark ap- 
plies to each of the systematic errors, namely that they 
cause an upward shift on the energy, of magnitude in- 
creasing with coverage (this behavior as a function of the 
coverage is an empirical result for the time step error, and 
an expected feature of all other sources of bias) . In par- 
ticular, the leading finite-size correction to the energy is 
given by the 2D Fermi energy contribution, which is pro- 
portional to the areal density and greater for the A^ = 26 
than for the infinite system; the population control bias 
PP] and the fixed-node error are both positive, and both 
vanish in the — > limit (where the trial wave function 
becomes exact). 

Thus, all of these errors combined are not expected to 
affect the main conclusion of the study, which is the pres- 
ence of a local minimum in the e{9) curve on Na and K 
sustrates. If anything, it is possible that one such mini- 
mum may exist for Li as well, as suggested by the results 
shown with open symbols in Figure [l] These energies 
are obtained by applying corrections to the simulation 
results, by estimating the magnitude of all systematic er- 
rors; in particular, the generally unknown magnitude of 
the fixed-node error is assessed by looking at the differ- 
ence between fixed-node and transient estimate [TU] for 
the purely 2D case. The correction is not large, and sug- 
gests that a liquid phase may exist on Li as well, with 
an equilibrium coverage of approximately 0.015 A^^ (a 
similar correction would not change quantitatively the 
physical conclusions for all other substrates). 

An interesting question is at what temperature the self- 
bound ^He film should turn superfluid. The methodology 
that we have adopted in this work does not allow us to 
explore this issue directly by simulation, due to the sign 
problem. An order of magnitude estimate might be ob- 
tained by computing the binding energy of a pair of "^He 
atoms in the presence of the substrate, on the assumption 
that superfluidity of a Fermi system should be underlain 
by the formation of bound pairs of atoms of opposite 
spins [3T1 132] ■ In Rcf. [28 , the claim is made that the 
binding energy, which is 20 /iK in 2D, could be as high 
as ~ 10 mK in the vicinity of a weak substrate. Thus, an 
order of magnitude estimate of the superfluid transition 
temperature is of the order of 1 mK. 

Summarizing, we have carried out first principle QMC 
simulations yielding robust numerical evidence to the ef- 
fect that ^He will form a stable, quasi-2D liquid phase at 
T=0 on Na, K, Rb and Cs substrates, and possibly on Li 



as well. The formation of the self-bound liquid phase is a 
consequence of the large atomic zero-point motion in the 
direction perpendicular to the substrate, which has the 
effect of softening the repulsive part of the helium inter- 
atomic potential at short distances. Based on this effect 
alone, no thin film is predicted to form on stronger sub- 
strates, such as Mg and graphite. The predictions made 
in this work appear to be experimentally testable, given 
the wealth of investigative work carried out over the past 
two decades, aimed at characterizing the physics of "^He 
films adsorbed on alkali substrates j3^Hi5] . 

The authors wish to acknowledge useful discussions 
with Milton W. Cole. 
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